Impurity band induced by point defects in graphene 
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It is pointed out that point defects on graphene are strongly correlated and can not be treated as 
independent scatters. In particular, for large on-site defect potential, it is shown that defects induce 
an impurity band with density of state characterized by the Wigner semi-circle law. We find that 
the impurity band enhances conductivity to the order of Ae^ /h and explains the absence of strong 
localization. Furthermore, the impurity band supports ferromagnetism with the induced magnetic 
moment approaching per defect in the limit of infinite quasi-particle lifetime. 
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Recent experimental realization of single-layer 
graphene [l|, Qhas raised much interest in studying 2D 
Dirac fermions in the context of condensed matter 
physics. One of graphene's peculiar properties is the 
anomalous electronic properties associated with defects 
and disorders. The notable example is the observed fer- 
romagnetic state induced by bombarding graphene with 
protons Another example is the finite conductivity 
at the Dirac point [4], indicating the absence of strong 
localization in graphene. These observations appear 
to deviate from what is expected for ideal and clean 
graphene. Since real graphene must involve disorders, it 
indicates that disorders may play an important role in 
these phenomena. 

Theoretically, it is known that the Dirac equation 
supports localized solutions in the presence of disor- 
ders. For instance, edge states would appear and are 
localized near the zig-zag edge[^. These localized so- 
lutions thus provide support for the observed mag- 
netism in carbon nanoribbons when the Coulomb inter- 
action is included[f|]. In general, magnetism due to de- 
fects in graphene gets widely support either by calcu- 
lations based on tight-binding models or first-principle 
calculations 0, H, [il. However, results obtained are vary- 
ing and appear to be sensitive to sizes and boundary 
conditions of the system ll| . It is not firmly established 
that ferromagnetism can exist in grapheneQ. On the 
other hand, theoretical work on the transport property 
of graphene also shows varying results on conductivity 
while experimental data consistently indicates that the 
conductivity is larger than the value Ae^ /nh obtained by 
many calculations It is therefore important to have a 
more reliable method for analyzing disordered graphene. 

Experimentally, there are many possible forms of de- 
fects in graphene^. [l3| . However, large defects with truly 
localized states interact via RKKY interaction and tend 
to induce antiferromagnetism[lo| . Thus these defects do 
not contribute the observed ferromagnetism. Theoreti- 
cally, it is known that the electronic state near a point 
defect or a disk-like defect on graphene is semi-localized 
with amplitude decaying as 1/r when the distance r of 



the electron to the defect is large IJ, |l5| . The semi- 
localization implies that disk-like defects interact directly 
and thus may be the source for ferromagnetism. Further- 
more, it explains the dependence on sizes and boundary 
conditions for numerical results. In particular, it implies 
that disk-like defects interact strongly and can not be 
treated perturbatively. Hence it calls for an appropriate 
method to take the semi-localization into consideration. 

In this Letter, we show that the semi-localization na- 
ture of defect states enables an impurity band form near 
zero energy. The existence of impurity band explains 
both the absence of localization and possible existence 
of ferromagnetism in graphene. Specifically, we shall in- 
vestigate an infinite graphene with randomly distributed 
point defects, characterized by an on-site potential u. In 
the large u limit [l^, it is shown that the electronic spec- 
trum due to randomly distributed point defects can be 
mapped into that of a random matrix. As a result, an 
impurity band with density of state characterized by the 
Wigner semi-circle law|17| forms. This impurity band 
appears to be observed in previous numerical results 
Due to the appearance of anomalous density of states 
near the Dirac point, we show that the conductivity 
is enhanced beyond Ae'^/nh . By further including the 
Coulomb interaction, it is shown that the impurity band 
supports ferromagnetism with induced magnetic moment 
depending on quasi-particle lifetime and defect density. 

We start by considering one point defect in graphene. 
Let the Hamiltonian for electrons in the tt band of an 
infinite graphene be Hq and the defect is located at r = 0. 
The wavefunction ip for an electron then satisfies 



(1) 



The solution ip can be found by resorting to the Huygens' 
principle which implies the existence of a particular so- 
lution propagating outwards from the defect to infinity. 
Obviously, this solution represents a state localized near 
the defect and is simply the Green's function, G(r, f ,E), 
which satisfies {E — Ho)G{'r, r' , E) = (5f.r'. Since both 
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and E are real, it is suffice to take a particular solution 

i^p{r)^ ARe[G{r,Q,E)], (2) 

where A is a, normalization constant. Clearly, because 
u6f^Qip{r) — u5f?o'0(O); to satisfy Eq.(IT]), one requires 
ASpQ — u5p_Qi}}p{r) — Au5pflRe[G{Q,Q,E)]. Hence the 
condition for the existence of the localized state at E^ is 
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- = Re [G(0,0,So)], 



(3) 



The above condition can be also reached by summing 
the Dyson series in which Gu{r,'r',E) = G{r,r',E) + 
G{r,0,E)G{Oy,E)/[l/u-G{0,0,E)] with G„ being the 
Green's function in the presence of the defect. Clearly, 
Eq obtained in Eq.Q specifies the pole in G„. In the 
limit u —f oo, only when Re[G{0,0, Eq)] — supports 
a localized state. For general u, the intersection of the 
curve Re [G(0, 0, E)] with 1/u gives the energy Eo of the 
localized state. Direct numerical calculations show that 
near E ~ 0, Re [G(0, 0, E)] w —^E. Hence for finite and 
large m, Eq ^ 1/u and in the limit u —^ oo, Eq approaches 
0. Here 7 characterizes the slope of the Green's function 
near E = 0. To the leading order of inverse quasi-particle 
lifetime, S, one finds 
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(4) 



(2^)2<52 + |^^|2' 

where ±|i5fe| are the energy spectra of the tt band. It 
shows that 7 goes to infinity as S approach zero. 

Note that ijjp is a particular solution with the bound- 
ary condition: ■0 — > as r ^ 00. Specifically, i^p{r) goes 
as 1 /r when r — > 00 IJ, |15| . Since 1 /r is not square- 
integrable, V'p is known as a semi-localized state. The 
semi- localization implies the solution to Eq.([T]) depends 
on the boundary condition. In general, the localized state 
hybridizes with extended states at the same energy. The 
degree of hybridization is determined by the boundary 
condition. In the case when -0 — > 0^ as r — s- 00 with (^fe 
being the plane wave function with Ek = Eq, by appro- 
priate superposition of (j)k and Im[G{r,0, E)], one ob- 
tains ■(/'(r) = (j)k{r) + uG{r,0, Eo)tp{0), which reproduces 
the usual scattering solution obtained by the Lippmann- 
Schwinger equation. In this case, after expressing ip{0) 
by (pkiQ), it is clear that ipi^r) is still dominated by the 
semi-localized part, G{f,0, Eq). 

As indicated in the above, since G(r, 0, Eq) dominates 
in the weight even if it can hybridize with extended 
states, it is therefore sufficient to consider ijjp. In partic- 
ular, it implies that defects would couple strongly when 
many defects are in present. The strong coupling im- 
plies that point defects can not be treated perturbatively. 
Clearly, the above construction can be generalized to the 
case with Nj impurities located at fi with i = 1,2,3, ■■■ , 
and Nj. In this case, we have the obvious solution 



Ni 



(5) 



As a generalization of Eq. ([3]) , the energy of semi- localized 
states satisfy 
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~9Ni1 ~9Ni2 ~9Ni3 ■ ■ i/u — gNiNi 

where gij = [ReG{fi,fj, E)] — gji. For low density of 
defects, one expects that energies of semi-localized states 
are near 0. Since Re [G(0, 0, E)] w -jE for S - 0, in 
the limit of u — s- 00, Eq.Q reduces 
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where gf^ = Re[G{'ri,rj,0)]. For randomly distributed 
positions of defects, Eq.([7]) implies that except for the 
scaling factor 7, energies of the electronic states arc ex- 
actly the eigenvalues of a symmetric random matrix. In 
the limit iV/ ^ 00, it is known that the density of jE 
follows the Wigner semi-circle lawpLSj. Converting the 
Wigner semi-circle law to the density of electronic states, 
D{E), we find 



D{E) = ^V^^'-l'E^ 



(8) 



where nj is the density of defects and = Ni{i9ij)^) 
with (• • •) denoting the average over positions of defects. 
Thus the electronic states due to defects form an impu- 
rity band with width w being determined by A via the 
relation w — 4A/7. We note in passing that Eq.® can 
be also obtained by summing the Dyson series for gij. 
However, the existence of the impurity band is a non- 
perturbative result which can not be obtained by sum- 
ming each averaged term (averaged over the impurity 
potential) in the Dyson series. 

For randomly distributed defects on graphene with N 
lattice points, the probability for a defect being located 
at r, is 1/N. Hence < (gO-)^ >= E^J(4)V^^ which 
shows that < (5^^)^ > is nothing but the Fourier trans- 
formation of {gij)'^ at fc = 0. We find A^ — jnj/2 and 
hence w = y^Sni/^, which vanishes when 5 goes to zero. 
Note that because g° does not vanish only when i and j 
belong to different sublattices, the impurity band results 
from the bipartite nature of graphene. 

We first address the issue on possibility of magnetism 
induced by the impurity band. We note that the impurity 
band is a flat-band with hopping connecting every pair 
of defects. According to the mechanism for flat-band 
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FIG. 1: Locations of defects that favor the spin-triplet state 
when two electrons filled in semi-localized state of two defects. 
Here one defect is fixed at the center, marked by X, while 
the other defect is denoted by a solid dot only when the spin- 
triplet is favored. Here sizes of the dots represent the strength 
of the triplet state. 



ferromagnetism'l9], the impurity band supports ferro- 
magnetism for the case when the system can be modeled 
by the Hubbard model. However, since screening is con- 
sidered to be less severe in two dimension, it calls for 
using more realistic Coulomb potentials. To investigate 
the magnetic property due to realistic Coulomb poten- 
tial, we calculate the exchange energy for the impurity 
band due to unscreened long-range Coulomb interaction 
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We start by considering an undoped graphene with two 
defects, in which semi-localized states are filled with two 
electrons. In this case, solutions to Eq.([7]) split into isl- 
and E^. with the splitting, — = 2\g'^{fi — r*,)|. 
Hence the splitting oscillates in the same way as does 
and vanishes when i and j belong to the same sublat- 
tice. Obviously, two electrons can be filled in E± sep- 
arately or filled in the same states with the spin state 
being singlet. If the Hund's rule dominates, two elec- 
trons fill in E± separately with their spins being in the 
triplet state. Based on solutions to Eq.® and the cor- 
responding wavefunctions, both Coulomb and exchange 
integrals are calculated so that the preferred spin state 
can be found. Fig. 1 shows our numerical calculations on 
locations of defects that favor the spin-triplet state for 
u = oo. Here one defect is fixed at the center, marked 
by X, while the other defect is shown by a solid dot only 
when the spin-triplet is favored. It is seen that for large 
distances when the splitting is small, spin triplet state 
is always favored; while the spin singlet state is favored 
only for some of lattice points at short distances when 
the splitting is large and exchange energy gain does not 
win over the splitting. 



To extend the above analysis to finite density of de- 
fects, we first note that the normalization of the wave- 
function in Eq.® implies that (X)f V'p(^)) = 1- If we 
set i? ~ 0, we find {A^^A'^) — Sij/Njj. By expanding 
= J^Ej '^E9{n,rj, E)Cl;^ and approximating E to 
in g{ri, rj,E), it is clear that (Cj^Cia) is determined by 
(C^^Cbct)- We find that the competition between ferro- 
magnetic, anti-ferromagnetic, and nonmagnetic states is 
determined by the exchange energy 
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(10) 

where — N^/Ni are fractions of electrons in the state 
with the z component of spin being a. In the ferromag- 
netic state, n| ^ n^. Therefore, the Fermi energies for 
different spin states are also different. Let E„ be the 
corresponding Fermi energy in the spin state a. The to- 
tal energy in the impurity band for each spin state is 
f-w/2 dEED{E). For an undoped graphene, Ei — ~E-\. 

We find that — ?i| = 2{y^Jl — y'^ + sin^"'^ y)/7r with 
y — 'jE'i /2A, while the change of the total energy in the 
impurity band per side is 



Afc 



2nTW ' 



Btt 



1 



(1 



y 



2^3/2 



(11) 



Using Eg. pi))) , the dependence of the exchange energy 
per site on — can be extracted and we find 
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exchange 



{2+V3)nj [y^/l^ + sin^^ ; 



(12) 

where a — 2. 46 A for graphene. In deriving Eq. ()12[) . we 
have expressed g^^ and l/|ri — r^l in the momentum space 
and approximated functions that are smooth in q by the 
corresponding values at Dirac points. It is clear that min- 
imizing Ak-\- AE exchange with rcspcct to y determines the 
induced magnetic moment. Note that the value j/o where 
the minimum occurs is determined by the intersection 
of y\Jl — y^ + sin^^ y and sy with s = Q .Q12w{eV) / nj . 
Fig. 2(a) shows a typical dependence of magnetic moment 
per defect (mfiB) on the width (w) of the impurity band 
for fixed defect density. It is clear that increasing the 
bandwidth would reduce the induced magnetic moment, 
which is consistent with the trend shown in Fig. 1. In the 
limit of infinite quasi-particle lifetime, because 7 goes to 
infinity and w — y^Sni/j, s approaches zero. Hence yo 
approaches one so that the induced magnetic moment 
approaches 1/is per defect, in consistent with the first 
principle calculation _8|. On the contrary, when temper- 
ature increases, S increases and 7 decreases, hence the 
induced magnetic moment decreases. For fixed quasi- 
particle lifetime, because Afc and AEexchange have dif- 
ferent nj dependence: Afc '-^ ri'J'^ and AEf-^change nj, 
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FIG. 2: (a) Typical dependence of magnetic moment (mfiB) 
per defect on the width (w) of the impurity band for n/ = 
O.OOOl.(b) Typical dependence of m on the defect density nj 
for fixed quasi-particle lifetime. Here 5 = O.OOli. Clearly, for 
fixed quasi-particle lifetime, a threshold in defect density is 
necessary to stabilize the ferromagnetic state. 

their competition leads to a threshold in the defect den- 
sity to induce ferromagnetism, as shown in Fig. 2(b). 

We now address the issue of the transport property by 
explicitly calculating the conductivity using the Kubo 
formula[20j. Here the total current is give by J = 
iet/TiJ^j s ^^j+Sa^j'^ with t being the hopping ampli- 
tude and S connecting nearest neighbors. In zero tem- 
perature, we find the dc conductivity is given by 

a = 2TThN {J%E')D{E)D{E')\j,^^ j,,^^ , (13) 

where we have expressed Cia in terms of Ceo- By ex- 
pressing Jqq in the momentum space, we obtain 



where / = ^ ^ sirv\kj2)\Eu\^ / {5^ + \Et,\')\ For 
small 5, because / ~ l/(27rt^(5^), we find a is larger 
than the value Ae'^/h by the factor w"^ / AttS"^ . Since typi- 
cal mean free path for graphene is about 30 nanometers, 
m ~ 10-5-10-6. Therefore, for 5/t ~ 10-^, we find that 
the enhance factor is from the order 0(1), in agreement 
with experimental observations [l^ . 

To summarize, we show that point defects on graphene 
are strongly correlated so that an impurity band near 
zero energy forms with the density of state characterized 
by the Wigner semi-circle law. It is shown that the im- 
purity band supports ferromagnetism with induced mag- 
netic moment depending on quasi-particle lifetime and 
defect density. Furthermore, we find that the induced 



magnetic moment approaches IfiB per defect in the limit 
of infinite quasi-particle lifetime and a threshold in defect 
density is required to stabilize ferromagnetism at fixed 
temperatures. Finally, the impurity band enhances the 
conductivity at the Dirac point to the order of Ae^/h, in 
consistent with experimental observations. 
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